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SUMMARY 


The  forces  and  overturning  moments  exerted  by  waves  on  large 
vertical  circular  cylinders  have  been  measured  in  the  laboratory.  Two  rigid 
cylinders,  12  in.  and  26.5  in.  in  diameter,  extending  from  the  bottom  of  a 
wave  flume  through  the  water  surface,  were  tested  in  varying  depths  of  water, 
for  a range  of  wave  periods  and  wave  heights  up  to  the  point  of  breaking.  A 
digital  computer  was  used  for  the  acquisition,  processing,  plotting  and  storage 
of  the  experimental  data. 

In  addition  to  the  experimental  work,  a design  method  is  presented 
which  allows  the  wave  loads  on  large  circular  cylinders  to  be  estimated  by 
means  of  a simple  desk  calculation.  The  experimental  data  shows  that  this 
simple  method  of  calculation,  based  on  the  linear  diffraction  theory  of 
MacCamy  and  Fuchs,  is  accurate  over  a wide  range  of  wave  conditions  and 
cylinder  sizes. 


RESUME 


Les  forces  et  les  moments  de  renversement  exerces  par  les  houles 
sur  des  cylindres  larges,  circulates  et  verticaux  ont  ete  mesures  dans  le  labo- 
rs to  ire. 


Deux  cylindres  rigides  de  12  et  26.5  pouces  de  diam&tre,  fixes  au 
fond,  s’etendant  au-dela  de  la  surface  libre,  ont  ete  etudies  dans  des  diffe- 
rentes  profondeurs  d’eau.  On  a fait  varier  les  periodes  et  les  amplitudes  des 
houles  dans  une  large  gamme  jusqu’aux  limites  du  deferlement. 


Un  ordinateur  digital  a ete  employe  pour  l’acquisition,  l’analyse, 
le  tracage  et  le  stockage  de  donnees  experimentaux. 


En  plus  de  l’aspect  experimental,  on  presente  egalement  une 
methode  devaluation  qui  permet  de  determiner  les  forces  des  houles  sur  des 
cylindres  larges  et  circulates,  par  un  simple  calcul.  Les  essais  mettent  en 
evidence  que  cette  simple  methode  de  calcul  qui  est  basee  sur  la  theorie 
1 innate  de  diffraction  de  MacCamy  et  de  Fuchs  se  prete  3 un  large  eventail 
des  param&tres  de  houle  et  egalement  a des  differentes  dimensions  de  cy- 
lindre. 
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WAVE  LOADS  ON  LARGE  CIRCULAR  CYLINDERS: 
A DESIGN  METHOD 


1.0  INTRODUCTION 

When  considering  wave  forces  on  rigid  circular  cylinders,  it  is  necessary  to  be  aware  of  the 
differences  between  the  wave  forces  on  large  and  small  diameter  cylinders.  If  the  diameter  of  a vertical 
cylinder  is  small  relative  to  the  wavelength,  the  wave  force  on  the  cylinder  is  normally  calculated  using 
the  Morison  equation  (Morison  et  al.1 3 ) which  consists  of  the  sum  of  the  inertial  force  and  the  viscous 
drag  force.  The  Morison  equation  simplifies  the  general  problem  of  wave  structure  interaction  by  as- 
suming that  the  cylinder  does  not  disturb  the  incident  wave  in  any  way.  But  as  the  diameter  becomes 
large  relative  to  the  incident  wavelength,  the  large  diameter  cylinder  causes  reflection  and  diffraction 
of  the  incident  waves,  the  Morison  equation  does  not  apply  and  a diffraction  theory  must  be  used. 
Viscous  drag  forces  are  not  important  for  wave  heights  which  are  small  relative  to  the  cylinder  diameter, 
which  is  the  usual  case  for  large  diameter  cylinders  and  thus  inertial  forces  are  always  predominant. 

MacCamy  and  Fuchs1  1 developed  a diffraction  theory  for  calculating  wave  loads  on  a vertical 
circular  cylinder  extending  from  above  the  water  surface  to  the  mudline;  however,  there  are  few  exper- 
imental results  to  confirm  their  theory.  Jen8  has  presented  data  agreeing  with  the  theory  but  it  is  only 
for  a very  limited  range  of  cylinder  diameter  on  wavelength  D/L  up  to  0.12.  The  experimental  results 
of  Lebreton  and  Cormault1 0 are  only  for  D/L  equal  to  0.09  and  do  not  compare  favourably  with  the 
diffraction  theory.  Data  has  been  published  by  Watanabe  and  Horikawa1 7 for  D/L  up  to  0.4,  but  there 
are  significant  differences  between  the  experimental  results  and  the  diffraction  theory.  Experiments 
by  Chakrabarti  and  Tam3  for  D/L  up  to  0.54  and  one  test  at  D/L  of  1.3,  show  agreement  with  MacCamy 
and  Fuchs’  theory;  however,  as  with  the  other  data  there  is  very  little  information  concerning  the  appli- 
cability of  the  theory  for  large  wave  steepnesses. 

Chakrabarti1  obtained  an  expression  for  the  wave  force  on  a cylinder  due  to  Stokes’  fifth- 
order  wave  theory,  but  the  kinematic  free  surface  boundary  condition  in  the  vicinity  of  the  cylinder 
was  not  satisfied.  Chakrabarti2  later  compared  this  approximate  theory  limited  to  the  second  order, 
to  a second-order  theory  obtained  using  a perturbation  method  by  Yamaguchi  and  Tsuchiya1 8 . He 
found  that  although  the  experimental  data  obtained  by  Yamaguchi  and  Tsuchiya1 8 correlated  well 
with  their  second-order  solution,  his  theory  and  experimental  data  showed  that  their  solution  over- 
estimated the  forces  considerably.  More  recently,  Raman  and  Venkatanarasaiah 1 5 have  also  used  a 
perturbation  method  to  develop  a second-order  diffraction  theory.  Unfortunately,  the  experimental 
results  given  to  confirm  the  theory  were  obtained  using  a relatively  small  diameter  cylinder  for  which 
previous  authors  have  assumed  that  viscous  drag  forces  are  significant.  The  results  would  have  been 
more  convincing  had  the  tests  been  conducted  on  a larger  cylinder  or  had  an  attempt  been  made  to 
separate  the  effects  of  viscosity  and  nonlinearity  of  the  waves. 

Wave  loads  on  large  floating  and  submerged  bodies  of  different  shapes  have  been  determined 
using  numerical  techniques  by  Garrison  and  Chow6 , Hogben  and  Standing7 , Lebreton  and  Cormault1  °, 
Garrett5 , Kokkinowrachos  and  Wilckens9  and  Van  Oortmerssen1 6 . Although  these  methods  can  be 
used  to  calculate  the  wave  loads  on  circular  cylinders,  the  computer  programs  are  time  consuming  to 
develop  and  expensive  to  use.  This  report  presents  a design  method  based  on  the  linear  diffraction 
theory  of  MacCamy  and  Fuchs1 1 to  estimate  wave  forces  and  moments  on  large  vertical  circular 
cylinders.  The  method  has  been  simplified  so  that  only  three  graphs  need  be  used  to  obtain  the 
complete  solution.  Force  and  moment  measurements  on  circular  cylinders  show  that  the  theory  gives 
a satisfactory  solution  over  a large  range  of  wave  conditions  and  relative  cylinder  sizes. 

2.0  THEORETICAL  METHOD 

Is  is  assumed  that  the  fluid  is  inviscid  and  homogeneous  and  that  the  fluid  motion  is  irrota- 
tional.  The  heights  of  the  waves  are  small  relative  to  the  wavelengths  so  that  terms  involving  wave 
steepness  to  second  or  higher  order  powers  may  be  neglected. 
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The  co-ordinates  are  chosen  such  that  x is  positive  in  the  direction  of  wave  propagation  and 
y is  positive  in  the  upward  direction.  The  origin  of  these  co-ordinates  is  at  still  water  level  where  the 
water  surface  elevation  plotted  against  x has  a maximum  negative  slope  (Fig.  1). 

According  to  small  amplitude  wave  theory,  the  velocity  potential  of  the  waves  incident  on 
the  rigid  cylinder  is 


<t>\  = - 


gH 

2o 


COSk<y  1*1  gitkx  - a,) 

cosh  kd 


(1) 


where  g is  the  acceleration  due  to  gravity,  H is  the  wave  height,  a is  the  wave  frequency,  a = 2n/T, 

T is  the  wave  period,  k is  the  wave  number,  k = 2ir/L,  L is  the  wavelength,  d is  the  water  depth,  t is  the 
time.  The  equation  for  plane  incident  waves  can  be  expressed  as  an  infinite  series  using  the  polar  co- 
ordinates r and  0: 


gH  cosh  k(y  + d) 
2o  cosh  kd 


2 em  im  cos  m0  Jm  (kr)  e_iot 
m * 0 


(2) 


where  Jm  (kr)  are  Bessel  functions  of  the  first  kind  of  orders  0, 1,  2 . . .m,  e0  = 1 and  em  = 2 for 

m > 1. 

A cylindrical  wave  is  reflected  from  the  cylinder  and  may  be  described  by  the  velocity 
potential  0R : 


0R  ~ 


gH  cosh  k(y  + d) 
2o  cosh  kd 


2 An,  cos  m0  [Jm  (kr)  + i Ym  (kr)]  e-iot 
m = 0 


(3) 


where  Am  are  constants  and  Ym  (kr)  are  Bessel  functions  of  the  second  kind. 

The  total  potential  <t>  is  the  sum  of  the  incident  and  reflected  potentials: 


gH  cosh  k(y  + d) 

2o  cosh  kd  m = o 


2 [em  im  cos  m0  Jm  (kr)  + Am  cosmfl  H^  (kr)|  e‘io1 


(4) 


where  H^  ’(kr)  = Jm  (kr)  + i Ym  (kr)  is  the  Hankel  function  of  the  first  kind.  It  should  be  noted  that 
MacCamy  and  Fuchs1 1 in  their  equation  for  the  total  potential,  have  incorrectly  used  the  Hankel 
function  of  the  second  kind  which  in  combination  with  e_io1  defines  a cylindrical  wave  moving  towards 
the  cylinder.  The  constants  Am  are  evaluated  by  using  the  boundary  condition  that  the  fluid  velocity 
normal  to  the  cylinder  is  zero,  that  is, 


— = 0 at  r = D/2  = a (5) 

dr 

where  D and  a are  the  diameter  and  radius  of  the  cylinder  respectively.  Thus, 


■C  (ka) 
Hm* ) (ka) 


Am  “ em  *m 


(6) 
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for  m = 0,  1,  2,  3 . . where  J|„  (ka)  and  H^1  * (ka)  are  the  derivatives  of  Jm  (kr)  and  ’ (kr)  at  r = a. 
Therefore,  the  total  potential  is 


gH  cosh  k(y  + d) 
2o  cosh  kd 


I emP[Jm(kr) 

m = 0 


J-(ka)  m 

H*‘  ’(kr)l  cos  mO  e-'°{ 

H^>(ka) 


(7) 


The  water  surface  elevation  r?  and  the  dynamic  pressure  at  the  surface  of  the  cylinder  are  calculated 
using  the  total  velocity  potential  and  the  linear  Bernoulli  equation: 


p = -pgy 


(8) 


where  p is  the  fluid  pressure  and  p is  the  mass  density  of  the  fluid.  At  the  water  surface  p 0 and 
y = rj,  so  from  Equation  (8), 


V = 


H 

2 


£ emim  + 1 [Jm  (kr) 

m = 0 


J,'n  (ka) 

Hm"(ka) 


H * ( kr ) ) cos  m 0 e_lot 


Similarly,  the  dynamic  pressure  is 


(9) 


P 


pgH  cosh  k(y  + d) 
2 cosh  kd 


e im  + 1 

2 — [Jm  (ka)  H^n  (ka)  - (ka)  H^*  (ka)]  cosmO  e_iot 

m = 0 H^u(ka) 

(10) 


which  reduces  to 


pgH  cosh  k(y  + d)  " 

p — 2, 


i 


7rka  cosh  kd  m = o H'(1,(ka) 


cos  mO  e~,ot 


(ID 


The  horizontal  force  in  the  longitudinal  direction  per  unit  length  of  cylinder  is  calculated  by 
integrating  the  pressure  given  by  Equation  (11)  around  the  cylinder: 


2pgH  cosh  k(y  + d)  J “ 

fv  = — : — TT-:  / ^ 


€ im 

cm  1 


cos  mO  cos 0 dO  e~lot 


Trk  cosh  kd  o m » o H^u(ka) 
After  integration  this  becomes 


(12) 


2pgH  cosh  k(y  + d) 
k cosh  kd 


I ; -lot 


|J'i  (ka)  + i Y’,  (ka)]  i e 


(13) 


and  taking  the  real  part  only, 


2pgH  cosh  k(y  + d) 

f = A(ka)cos(ot-  a) 

k cosh  kd 


(14) 
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where 


and 


A(ka)  = |J',2(ka)  + Y',2(ka)r1/2 

/J'l  (ka)  \ 

\Y',  (ka )/ 


a = tan"  1 


The  total  horizontal  force  on  the  cylinder  is  obtained  by  integrating  Equation  (14)  from  the  bottom 
y = - d,  to  the  still  water  level  y = 0,  in  keeping  with  linear  theory: 


F(t)  = 


2pgH 


A(ka)  tanh  kd  cos  (at  - a) 


(15) 


The  overturning  moment  about  the  base  of  the  cylinder  is 


M(t)  = / f„  (y  + d)  dy 
-d 


Substituting  for  fx  from  Equation  (14)  and  carrying  out  the  integration  gives 


(16) 


M(t) 


2pgH 


A(ka)  |kd  tanh  kd  + sech  kd  - 1]  cos  (at  - a) 


(17) 


The  distance  yt.  measured  from  the  bottom  to  the  point  at  which  the  resultant  horizontal  force  acts  is 
obtained  by  dividing  Equation  (17)  by  Equation  (15): 


yc 


= — | kd  + csch  kd  - coth  kd] 
2tt 


(18) 


The  horizontal  force  per  unit  length  given  by  Equation  (14)  may  be  written  equivalently 
as  an  inertia  force: 


pjrD2  du 
4 dt 


(19) 


where  Cm  is  the  coefficient  of  mass  and  du/dt  is  the  horizontal  component  of  water  particle  accelera- 
tion. Therefore, 


Ptt3D2H  coshk(y  + d) 

Cm  — . . . cos  (kx  - ot) 


2T2 


sinh  kd 


(20) 


Equating  Equations  (14)  and  (20)  gives  an  expression  for  Cm  at  x = 0 at  the  centre  of  the  cylinder: 


4L2  cos (at  - a) 

Cm  = A(ka) 


m tt3D2 


cos  ot 


(21) 


A modified  coefficient  of  mass  C*  is  defined  as 


5- 


C, 


41/ 
IT3  D2 


A(ka) 


(22) 


such  that 


cos  (at  - a) 
cos  at 


(23) 


The  coefficient  of  mass  Cm  is  a function  of  time  and  varies  through  a wave  cycle  as  shown  in  Figure  2. 
It  is  constant  through  the  wave  cycle  only  when  a is  zero,  that  is,  for  D/L  approaching  zero  or  approx- 
imately 0.59.  When  a is  not  zero,  Cm  is  equal  to  C*  when 


at  = tan" 1 


( 


1 - cos  a 
sin  a 


) 


(24) 


The  total  force  in  the  longitudinal  direction  is  obtained  by  integrating  the  force  per  unit 
length  as  given  by  Equation  (20)  for  x = 0,  through  the  depth  of  the  water  to  still  water  level: 


F(t)  = Cn 


P7T2  D2  HL 
4T2 


cos  ot 


(25) 


Substituting  the  expression  for  Cm  given  by  Equation  (23)  into  Equation  (25)  gives 


F(t)  = C* 


f)Tt2  D2  HL 
m 4^2 


cos (at  - a) 


(26) 


The  maximum  force  occurs  when  at  = a: 


Fjiiax 


_ , pit2  D2  HL 


4T2 


(27) 


The  overturning  moment  about  the  base  of  the  cylinder  can  also  be  expressed  in  terms  of  C*,  as 


M(t)  = C, 


pgHLD2 


16 


|kd  tanh  kd  + sech  kd  - 1)  cos  (ot  - a) 


(28) 


The  maximum  overturning  moment  occurs  when  at  = a: 


M 


mux 


,*  PgHLD2 

m 1C 


| kd  tanh  kd  + sech  kd  - 1 1 


(29) 


Figures  3,  4 and  5 are  design  curves  based  on  the  above  equations.  Figure  3 gives  C*  by  the 
solution  of  Equation  (22),  so  that  the  force  Fmax  may  be  calculated  by  using  Equation  (27).  The 
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overturning  moment  about  the  base  of  the  cylinder  Mmax  given  by  Equation  (29),  can  be  obtained  in 
Figure  4.  The  phase  angle  a is  plotted  in  Figure  5 and  allows  the  determination  of  F(t)  and  M(t)  given 
by  Equations  (26)  and  (28)  respectively.  Moments  may  also  be  calculated  by  the  product  of  the  force 
and  lever  arm  (Fig.  6). 


3.0  EXPERIMENTAL  METHOD 


The  circular  cylinder  (Fig.  7)  used  in  the  experiments  was  12  in.  in  diameter  and  was  con- 
structed of  1/4  in.  thick  plexiglass.  It  was  supported  1/8  in.  above  the  wave  flume  bottom  by  a rigid 
3 in.  diameter  steel  tube  clamped  to  a steel  frame  above  the  flume.  A force  meter  was  contained  within 
the  cylinder  and  consisted  of  two  3/4  in.  diameter  stainless  steel  strain  rods  12  in.  long.  The  wave  force 
on  the  cylinder  was  transferred  to  the  top  strain  rod  through  a horizontal  steel  bar  and  to  the  bottom 
strain  rod  through  a circular  steel  base  plate  (Fig.  7).  Foil  strain  gauges  glued  on  the  strain  rods  were 
aligned  so  that  total  horizontal  forces  could  be  measured  in  the  direction  of  wave  propagation  and 
normal  to  the  wave  direction  to  give  longitudinal  and  transverse  forces  respectively.  The  lower  strain 
rod  is  visible  in  Figure  7,  but  the  upper  strain  rod  is  hidden  inside  the  3 in.  diameter  tube.  Also  inside 
the  tube  are  two  solid  bushings  which  fix  one  end  of  each  of  the  strain  rods.  Using  the  Wheatstone 
bridge  outputs  from  the  upper  and  lower  strain  rods,  it  was  possible  to  measure  total  forces  and  to 
calculate  the  corresponding  overturning  moments  about  the  base  of  the  cylinder.  A more  detailed 
description  of  the  wave  force  meter  and  its  calibration  is  given  by  Pratte  et  al. 1 4 . From  the  calibration 
curves  of  the  force  meter,  its  error  band  was  estimated  to  be  less  than  ±2%  over  the  range  of  forces 
measured.  The  natural  frequency  of  vibration  of  the  cylinder  in  the  maximum  depth  of  water  was 
approximately  11  Hz. 


Wave  profiles  in  the  flume  were  measured  by  the  noncontacting  capacitive  wave  transducer 
(Zwarts1  9 .2  0 ) shown  in  Figure  7.  It  was  suspended  above  the  water  surface,  midway  between  the 
cylinder  and  the  flume  wall  where  the  wave  heights  measured  were  approximately  equal  to  the  incident 
wave  heights  measured  without  the  cylinder  in  the  flume.  The  wave  flume  was  approximately  6 ft.  wide, 
4.5  ft.  deep  and  220  ft.  long. 


A limited  number  of  tests  were  conducted  using  a 26.5  in.  diameter  cylinder  to  obtain  force 
and  moment  data  for  larger  values  of  relative  diameter  D/L.  The  cylinder  was  suspended  from  a force 
meter  in  a flume  12  ft.  wide,  4.5  ft.  deep  and  162  ft.  long.  The  force  meter,  shown  in  Figure  8 with  a 
square-section  cylinder  attached  to  it,  consisted  of  three  aluminum  strain  members  3 in.  in  diameter. 
The  strains  in  these  members  caused  by  wave  loads  were  measured  using  semiconductor  strain  gauges 
forming  six  Wheatstone  bridges.  Three  bridges  measured  strain  due  to  bending  and  gave  outputs  pro- 
portional to  the  longitudinal,  transverse  and  vertical  forces.  Three  bridges  measured  strain  due  to 
torque  and  gave  outputs  proportional  to  the  moments  about  the  three  co-ordinate  axes.  A detailed 
description  of  this  force  meter  is  given  by  Funke4 . Wave  heights  were  measured  using  a capacitive 
rod  wave  gauge  at  the  location  of  the  model  but  withsqt  the  model  in  the  flume. 


Monochromatic  waves  were  generated  in  five  water  depths  ranging  from  9.7  in.  to  29.0  in. 
and  seven  wave  periods  were  tested  from  0.77  sec.  to  2.58  sec.  For  each  water  depth  and  period,  a 
number  of  waves  were  generated,  ranging  from  small  amplitude  waves  to  those  that  were  on  the  point 
of  breaking. 


A digital  computer  was  used  for  the  acquisition,  processing,  plotting  and  storage  of  data. 
Sampling  of  the  data  was  commenced  after  the  waves  passing  the  model  reached  steady  state  conditions. 
For  each  test,  the  wave  profile  and  the  corresponding  wave  forces  were  sampled  every  one  hundredth 
of  a second  for  a total  of  six  seconds.  Total  forces  and  overturning  moments  were  computed  and  auto- 
matically plotted  along  with  the  wave  profile  measured  at  the  cylinder.  The  experimental  results  were 
also  printed  and  the  data  was  stored  on  magnetic  tape. 
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4.0  EXPERIMENTAL  RESULTS 

The  number  of  waves  in  each  test  record  varied  depending  on  the  wave  period  because  of  the 
fixed  sampling  time  of  six  seconds.  Any  test  with  a variation  in  wave  height  or  maximum  force  meas- 
urement within  the  record  of  more  than  5%  was  not  used.  From  the  test  records,  the  forces  Fniax  and 
overturning  moments  Mmax  were  taken  as  average  absolute  maximum  values  of  positive  and  negative 
measurements.  The  difference  between  the  positive  values  and  the  averages  were  normally  less  than  6%. 
The  forces  and  corresponding  moments  measured  in  the  transverse  direction  were  negligible  and  are 
not  included  in  the  presentation  of  experimental  results. 

The  maximum  forces  and  moments  were  expressed  nondimensionally  as  functions  of  wave 
steepness  H/L,  relative  diameter  D/L  and  relative  depth  d/L.  This  form  of  data  presentation  is  used  to 
show  the  variation  of  wave  forces  and  moments  with  wave  steepness.  The  graphs  in  Figures  9 to  22 
show  the  comparison  between  the  experimental  results  and  the  theoretical  forces  and  moments  cal- 
culated by  Equations  (27)  and  (29).  In  Figures  9 to  18  for  which  the  relative  diameter  D/L  is  greater 
than  0.09,  there  is  close  agreement  between  the  theory  and  experimental  results.  However,  for  D/L 
less  than  0.09,  the  measured  data  differs  considerably  from  the  theory  particularly  for  very  small  D/L 
and  large  wave  steepness  H/L.  For  this  data,  the  absolute  values  of  positive  and  negative  wave  loads 
are  no  longer  approximately  equal.  They  differ  from  their  average  values  by  as  much  as  8%  for  the 
data  in  Figure  20  for  which  D/L  = 0.08  and  20%  for  the  data  in  Figure  22  for  which  D/L  = 0.057. 

In  addition,  for  large  wave  steepnesses  the  phase  angles  between  the  wave  load  records  and  the  wave 
profiles  do  not  correspond  to  the  theoretical  values  for  a given  in  Equation  (14).  The  phase  angle 
increases  with  increasing  wave  steepness,  which  is  typical  of  the  effect  of  an  increasing  drag  force. 

These  characteristics  of  the  experimental  results  are  much  more  pronounced  in  the  measurements  of 
wave  loads  on  square-section  cylinders  as  described  by  Mogridge  and  Jamieson1 2 . 

Although  it  appears  that  the  deviations  from  theory  described  above  are  due  primarily  to 
viscous  drag  forces,  similar  deviations  will  occur  due  to  nonlinearity  of  the  waves.  The  lowest  value 
of  d/L  measured  was  0.09,  for  which  experimental  data  is  presented  in  Figure  16  and  also  Figure  22. 

In  Figure  22,  D/L  is  less  than  0.09  indicating  that  there  is  a viscosity  effect  but  the  deviation  of  the 
experimental  data  from  the  theory  could  also  be  caused  by  the  nonlinearity  of  the  waves.  The  fact 
that  it  is  not  caused  by  nonlinearity  is  evident  in  Figure  16  where  the  experimental  data  agrees  with 
the  theory  at  the  same  value  of  relative  depth  d/L  = 0.09.  Mogridge  and  Jamieson1 2 give  test  results 
for  square-section  cylinders  in  relative  water  depths  to  d/L  = 0.06  and  specify  the  limit  of  d/L  = 0.09 
at  which  nonlinearity  of  the  waves  causes  deviation  of  results  from  the  linear  theory. 

The  test  results  were  also  plotted  on  the  theoretical  design  curves  developed  from  the  linear 
diffraction  theory.  Using  the  measured  absolute  maximum  forces  with  known  values  of  wave  height, 
length,  period  and  cylinder  size,  C£  was  calculated  for  varying  values  of  D/L  and  plotted  on  the  theo- 
retical curve  in  Figure  23.  Dimensionless  moments  were  also  calculated  using  the  known  values  of 
C*  and  were  plotted  on  the  theoretical  curve  in  Figure  24.  The  data  used  is  the  average  of  a number 
of  tests  conducted  with  a range  of  wave  steepnesses.  For  tests  with  D/L  less  than  0.09,  the  data 
averaged  is  from  tests  with  wave  steepnesses  of  less  than  0.01.  For  D/L  greater  than  0.09,  the  data  has 
been  averaged  for  all  tests  run  at  each  value  of  D/L  and  d/L. 

The  experimental  data  plotted  in  Figures  23  and  24  show  good  agreement  with  the  linear 
diffraction  theory  for  D/L  between  0.057  and  0.441  and  also  for  d/L  between  0.090  and  0.786.  The 
maximum  wave  steepness  for  which  there  is  reasonable  agreement  between  theory  and  experiment, 
depends  mainly  on  the  relative  diameter  D/L.  Figures  9 to  22  indicate  that  above  a value  of  D/L  of 
0.090,  the  theory  may  be  used  for  wave  steepnesses  up  to  0.076.  However,  for  decreasing  values  of 
D/L,  the  theory  can  only  be  used  for  decreasing  wave  steepnesses,  until  at  D/L  equal  to  0.057  the 
error  at  H/L  equal  to  0.018  is  approximately  20%. 

5.0  CONCLUSIONS 

1.  A simple  design  method  based  on  the  linear  diffraction  theory  of  MacCamy  and  Fuchs  has  been 

developed  to  determine  the  wave  forces  and  overturning  moments  on  large  circular  cylinders. 
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2.  Results  of  experiments  have  shown  that  the  theoretical  method  can  be  used  over  a wide 

range  of  conditions;  that  is,  for  D/L  between  0.09  and  0.44,  for  d/L  between  0.09  and  0.79 
and  for  wave  steepnesses  up  to  0.08.  For  D/L  less  than  0.09,  the  theory  can  be  used  with 
reasonable  accuracy  for  low  wave  steepnesses. 
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FIG.  2:  VARIATION  OF  Cm  THROUGH  A WAVE  CYCLE 


- 13  - 


X 

cn 

Q-  10 
*E  p 
o 8 


FIG.  4:  DIMENSIONLESS  OVERTURNING  MOMENT  AS  A FUNCTION  OF  d/L 
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FIG.  8:  SIX  COMPONENT  FORCE  METER 
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FIG.  10:  MAXIMUM  FORCES  AND  MOMENTS 
(D/L  = 0.325,  d/L  = 0.786) 
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FIG.  11:  MAXIMUM  FORCES  AND  MOMENTS 
(D/L  = 0.239,  d/L  = 0.578) 
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FIG.  13:  MAXIMUM  FORCES  AND  MOMENTS 
(D/L  = 0.134,  d/L  = 0.213) 
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FIG.  19:  MAXIMUM  FORCES  AND  MOMENTS 
(D/L  « 0.086,  d/L  « 0.136) 
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FIG.  22:  MAXIMUM  FORCES  AND  MOMENTS 
(D/L  - 0.057,  d/L  « 0.090) 
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FIG.  24:  EXPERIMENTAL  RESULTS  FOR  OVERTURNING  MOMENTS  AS  A FUNCTION  OF  d/L 
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